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Pharmacokinetics for Regulatory Risk Analysis: 
The Case Of Trichloroethylene 


KF.NNETH T. BOGES 

Environmental Silences Division. Lawrence Livermore Sattoruti Laboratory. 
Livermore, California 945?0 
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Ph.vsialogictilh based pharmacokinetic ( PBPKI models describing the uptake, metabolism, 
and cccretion nf volatile organic compounds t VOTG are now proposed tor use in regulatory 
heallh-r.sk assessment. A steady-state analysis of one such model is shown to provide simple, 
convenient predicted relationships between an applied dose and the corresponding toxicologi¬ 
cal E effective, metabolized dose for certain VOCs like trichloroethylene I TCE). A version of 
this PBPK model was tit to data on human metabolism of TCE to urinary metabolites in chroni¬ 
cal!) exposed workers. > icidtne a direct estimate of PBPK parameters governing human capacity 
to metabolize TCP. Il is shown that this estimate is consistent with others based on experimental 
studies of TCE metabolism in humans exposed to TCE by inhalation for short periods. These 
results are applied to human cancer-risk assessment using rodent bioassay data on TCE-induced 
lumon^enesis. IdM Academic Press. Inc. 


INTRODUCTION 


Physiologically based pharmacokinetic (PBPK) models describing the uptake, me¬ 
tabolism, and excretion of volatile organic compounds ( VOCs) are now proposed for 
use in regulatory risk assessment. Here, the PBPK approach is used to estimate the 
extent to which humans metabolize trichloroethylene (TCE). TCE is One of the halo- 
genated VOCs whose metabolism, as opposed to applied dose, is directly correlated 
with experimentally induced chronic toxicity (Buben and O Flaherty. 1985) and I For 
regulatory purposes, presumptively) carcinogenicity (U.S. EPA, 1985). For compar¬ 
ative purposes, similar estimates based on short-term experimental studies of TCE 
uptake and metabolism in humans are first sum marized. and two early PBPK models 
applied to TCE arc reviewed. A more recent PBPK model developed by Ramsey and 
Anderson (1984) is then described and this model’s predicted behavior at steady state 
is analyzed. The steady-state analysis focuses on a quantity of key concern in risk 
analysis for VOCs like TCE. namely, the fraction of a potentially available dose that is 
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actually metabolized, particularly under conditions of low-level multiroute exposure 
relevant to environmental regulation. 

A version of the Ramsey-Andersen PBPK model is then tit to data on human 
metabolism of TCE to urinary metabolites in chronically exposed workers. The re¬ 
sulting estimate of the human capacity to metabolize TCE is shown to be consistent 
with earlier estimates based on the short-term, experimental studies referred to above. 
Finally, the steady-state analysis of the Ramsey-Andersen PBPK model is applied to 
an illustrative cancer-risk assessment for TCE, based on the former estimate of hu¬ 
man metabolic capacity and on estimates of TCE’s carcinogenic potency derived 
from rodem-bioassay data. 

EXPERIMENTAL STUDIES OP ICE METABOLISM IN HUMANS 


T( 


Study 
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Humans metabolize TCE primarily to trichloroacetic acid (TCA) and (bound and 
free) trichloroethanol (TCEL). which are excreted in urine for the most pail, al¬ 
though some TCEL and TCA may be further metabolized (Soucek and Vlachova, 
I960: Bartonicek. 1962. Nomiyama and Nomiyama, 1971: Muller etui.. 1972; Fer¬ 
nandez ei at.. 1975 . 1977; Monster et a!.. 1976. 1979). The extent to which people 
metabolize TCE during and after controlled respiratory exposures has been investi¬ 
gated by measuring total urinary metabolites generated postexposure or by quantify¬ 
ing the total amount of unchanged TCE expired after exposures for which total TCE 
retention was estimated {Soucek and Vlachova. 1960: Ogata et id. 1971; Nomiyama 
and Nomiyama. 1971; Ertlc et a!.. 1972; Muller et at. 1972: Monster et a!.. 1976; 
Fernandez el at., 1975. 1977; Monster et at.. 1979). Results from such studies are 
summarized in Tabic 1 The reported values for average percentage of retained dose 
metabolized (PRDM) range from 81 to 92%. Considering only the most recent stud¬ 
ies—those by Monster et at. and Fernandez et at. involving a total of 18 males given 
respiratory exposures ranging from 280 to 1440 ppm-hr TCE at TCE concentrations 
ranging from 54 to 140 ppm—the reported values for average PR DM are remarkably 
consistent, ranging from 89 to 92% with a person-weighted average of 90.6 ± 1.9% 
(CV = 2.1%) and a corresponding individual PRDM range of 87 to 93%. These 18 
individual PRDM values are not significantly correlated with either TCE exposure 
concentration or time-integrated TCE exposure (in ppm-hr). 


PHYSIOLOGICALLY BASED PHARMACOKINETIC 
MODELS APPLIED TO TCE 

Two early applications of the PBPK approach to TCE pharmacokinetics in hu¬ 
mans are reviewed briefly, followed by a description of a more recent PBPK model 
that has been applied to TCE. 


Two Early PBPK Models Applied to TCE 

A PBPK model developed by Sato et a!. (1977) was used to interpret data gathered 
on excretion of TCE and metabolites after respiratory exposure of four human volun¬ 
teers to 100 ppm TCE for 4 hr. A first-order PBPK model was used to simulate con¬ 
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TABLE : 


I CE Metabolism in Humans aeter Experimental Respiratory Exposure 


Si udy 

TCE exposure 

Collection 
period g 
(days) 

Sex (No.) 

Average % 
retained dose 
metabolized * - PRDM 
(range) 

Soucek and 

Vfachova ( I960) 

93-1 5^ ppm-5 hr 

10-14 

M<3> 

F O) 

— 

Ogata er at (1971) 

170 pprri'3 hr 

4.2 

?(4> 

— 


170 ppm-7 hr 

4.2 

7(5) 

— 

Nomivama and 
Nomiyama (1971) 

310 ppm-l60min 

6 

M(5) 

F (5) 

80.8 

87.3 

Enle eta!. (1972) 

50 ppm-6 hr for 5 days 

7 

M (5) 

— 


100 ppm-6 hr lor 5 class 

7 

M (5) 

— 

Mtiller ci a!. ( 19~ 2) 

50 ppm-6 hr for 5 days 

18 

Ml 5) 

— 

Monster ft al. (1976) 

71.) ppm-4 hr 

0 

M (4) 

90 

(59-93) 


UO ppm-4 hr 

4 

M (4) 

92 

(91-93) 

Fernandez a al 
(1975. 197?) 

54 ppm-8 hr 

16 

M (2) 

92 

(91.7-92.3) 


97 ppm-K hr 

21 

M (3) 

91 

(90.5-93.0) 

Monster el al M979i 

■’3 ppm-4 hr for 5 da»s 

6 

Ml?) 

89 

188-91) 


J Period of routine urine-collection. including exposure period. !n the Monster et al. (1976) Study, routine 
urine collection stopped at 66 hr. and a final sample was taken on Day 9 postexposure. 

' Here calculated as 1 00 T minus the reported percentage of retained TCE that is expired as TCE subse¬ 
quent to exposure, for lho.se studies reporting such a value. Retained TCti is defined i n these studies as the 
difference between inspired and expired TCE concentration times the volume of air respired. The retained 
fraction, /i, of TCE in alveolar air was reported to be about 78% (Monster et of. 1976: assuming an alveolar 
ventilation rate of 5.9 lilers/min). 7 [ to 78% (Fernandez et at.. 1975). and 70% (Monster et al-. 1979). 


centrations in three compartments: vessel-rich ( richly perfused) tissues, muscle (and 
other poorly perfused) tissues, and fat (very poorly perfused) tissues. In this model, 
the compartments are assumed to be interconnected and initially at an equilibrium 
determined by tissue volumes and blood/air and tissue/blood partition coefficients. 
The latter coefficients were based on experimental values obtained using rat and hu¬ 
man tissues. Richly perfused tissues are assumed to be the source of both metabolic 
and respiratory excretion of absorbed TCE. In this model, intercompartment ex¬ 
change of TCE is governed solely by postulated intertissue diffusion. When experi¬ 
mental human data were fitted to this model, a value of 104 iiters/hr was obtained 
for metabolic clearance of TCE. The authors did not explain the basis for their as¬ 
sumption that humans exposed to TCE would have reached compartmental equilib¬ 
rium after only 4 hr. 

A different, more complete PBPK. model was used by Fernandez el al (1977) to 
simulate the uptake, metabolism, and excretion of TCE in experimentally exposed 
humans. This model contains the three compartments used by Sato et at. (1977), 
plus a liver compartment and a pulmonary compartment. The model provides for 
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Fig. ! . Pathways jrid rale .’nr.sums Tor 1(1. metabolism used in the PBPK model of Fernandez r; u! 
(1977) describing uptake, metabolism, and excretion of TO: in humans, TCE. trichloroethylene. CH. 
chloralhydratc, TCEL. trichioroeihanoU TCEL-G, TCEL-glucuromde. 


mass balance of TCE in the compartments, TCE entering and leaving the body 
through the lung, and TCE excreted through metabolism in the liver to metabolites 
excreted in urine. The uptake and excretion kinetics were assumed to be linear and 
governed by tissue volumes, tissue-specific blood volumes, blood perfusions, and tis¬ 
sue/gas partition coefficients, in the model. TCE metabolism was assumed to take 
place solely in the liver via the pathways shown in Fig. 1. The values for the linear 
kinetic parameters used by Fernandez et al. (1977) were determined from the avail¬ 
able literature on TCE metabolism in humans, primarily from Muller ct al. (1974). 
These values are listed in Fig. 1. In this model. TCE metabolism is assumed to be 
limited only by blood flow through the liver in such a way that the fraction. / t . of 
arterial blood flowing into liver that is metabohcally cleared of TCE is a constant 
equal to 0.86. The liver-perfusion value of 96.3 liters/hr used by Fernandez et al. 
(1977) for a reference 70-kg man yields a metabolic clearance rate of 96.3 J c = 83 
liters/hr. which is close to the metabolic clearance rate of 104 liters/hr estimated by 
Sato et al. (1977), The PBPK model of Fernandez etal (1977) accurately predicted 
the decay in alveolar-TCE concentration observed over 50 hr and the cumulative 
urinary excretion of TCEL and TCA observed over 200 hr after volunteers, with an 
average weight of about 70 kg, were exposed for 8 hr to 54 to 160 ppm TCE. 

The Ramsev-Andersen PBPK Model 

The National Academy of Sciences (NAS) recently considered the use of PBPK 
models to facilitate dose-route extrapolation when using inhalation toxicity data to 
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Fig. 2. Schematic diagram of physiologically based pharmacokinetic (PBPK) model for inhalation of 
volatile organic compounds. The model assumes that four ''well-stirred'' compartments or tissue groups 
collect inhaled compound at rates governed by air concentration ( t' ln ). air and blood Hows I Q I. blood 
concentrations (ill. compartment volumes (I I. tissue/blood partition coefficients (PI, and metabolic 
paramete rs (! . A' m I. 


set safe drinking water limits (N’AS. 1986). A range of issues was considered in this 
study, including illustrative applications of the PBPK approach to dose-route extrap¬ 
olation from rats to humans for noncarcinogenic toxicity associated with exposure 
to TCE and benzene. The pharmacokinetic model used was developed by Ramsey 
and Andersen (I 984) to describe the uptake, metabolism, and excretion of styrene in 
rats and humans. The structure of the model is shown in Fig. 2, and its parameter 
definitions are given in Table 2, This type of model has been applied to the study 
and prediction of animal and human pharmacokinetics for other VOCs, including 
benzene, methylene chloride, and perchloroethylenet NAS. 1986; Gargase/ al., 1986; 
Andersen et al.. 1987; U.S. EPA, 1986a; Reitz and Nolan, 1986; Hattis et al., 1986; 
Ward et al., 1988). The model consists of a set of differential equations that describe 
the rate of change of the amount of absorbed chemical present in each of four physio¬ 
logically realistic tissue compartments, which are assumed to be ideally well-mixed 
at any given time. Metabolism is presumed to occur solely in the liver through a 
saturable enzymatic process with Michaelis-Menten kinetics. 

According to the Ramsey-Andersen model, pulmonary uptake of a VOC occurs 
continuously such that alveolar concentration. C a . is in instantaneous equilibrium 
with arterial blood concentration. 77,, governed by the blood/air partition coefficient. 
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table: 

Compartment and Parameter Di i-ini iions for the Ramshv 

-Andersen PBPK Model 

Abbrev. 

Definition 

Unit 

C,„ 

Concentration in air inhaled 

mg/liier air 


Concentralion in alveolar air 

mg/liier air 

{ exh 

Measured concentration in expired breath 

mg/litcr air 

a 

Alveolar ventilation rate 

liters air/hr 

Of 

Cardiac Output 

liters hlood/hr 

P„ 

Rlood/air partition coefficient 

titers air/litcr blood 

H, 

Arterial blood concentration 

mg/liier blood 

B 

Venous blood concentration 

mg/liter blood 


Amount metabolized in liver 

mg 

Q 

Blood flow rate to compartment i 

liters blood/'hr 

1 

Volume of compartment i 

liters i kg) 

c 

( once nt ration m compart mens i 

mg/liter 

H 

Concentration in venous blood leaving corn part ment i 

mg/liter blood 

.1 

Amount in compartment i 

mg 


Tissue/blood partition coefficient for compartment > 

liters blood ./liter tissue t 


Maximum metabolic rale 

mg/hr 

K m Miehaelisconstant 11, I n , a , _i: 

Com pa n m e n mis u bse nptx 
( l l i\er<metabolizing tissue «roup> 

^ \ at tissue i\crv pi>oi!\ ncitusedi 

3 Richl> perfused tissues (brain. kidne> „ uscerai 

mg/'liter blood 

4 Poorl 

> perfused tissues 1 muscle, skin) 



P h , in accordance wiih the relation C’ a = B 3 j P h . Similarly, the concentrations, C<, ot 
chemical in each tissue compartment are presumed to be in instantaneous equilib¬ 
rium with the concentrations. in venous blood exiting the corresponding tis¬ 
sue, governed by the corresponding tissue/blood partition coefficients such that B, 
= C./P, The amount of chemical in any given tissue compartment is given by 
A, - C'il ). For notational convenience, the dependence of state variables (Cs, B' s, 
and ,-fs) on time t is suppressed. 

1 n any given interval. (It. a VOC delivered to the lung via respiratory retention and 
via returning venous blood is balanced in this model by the chcmicaljnass exiting 
the lung via exhalation and via arterial blood, such that T Q b B]dt - {Q a C a 

+ Q b Bjdt. or at C n - CJ - Q h ( B a - B). which, recalling that C\ = BJP b , yields 


Qi, Cn + QbB 

(a/ a>) + Qb 


in 


(Note that an experimentally measured concentration in exhaled breath is given by 
C'ch = F d C, n + (1 - F d ) C a , where the dilution factor F d approaches zero as the effi¬ 
ciency of measuring purely alveolar expired air approaches 100%.) Equation (1) spec¬ 
ifies B 3 (and thus C B and C„h) to be at each instant a flow-weighted average of C.„ 
and B. Similarly, the concentration B in venous blood returning from each compart¬ 
ment is presumed to be the instantaneous flow-weighted average 
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( 2 ) 


For the nonmctabolizing tissues, the amount of chemical entering the ;th compart¬ 
ment via arterial blood during any given interval is set equal to the amount gained 
by that compartment plus the amount leaving in venous blood. The chemical con¬ 
centration in venous blood leaving each of these compartments is therefore defined 
by 


' = 2.3.4. (3) 

where dot notation is used here, and below, to represent differentiation with respect 
to time (i.e,. B, = dBJ dt. etc.). The amount of chemical metabolized in liver is given 
by the Michaclis-Mcnten relation 


- * max B \ 


(4) 


in which K m is defined as the concentration in venous blood from liver (or, alterna¬ 
tively, K m P i is the chemical concentration in liver) at which the liver’s metabolic 
velocity .-] m is onc-hulf its maximum value. I . The state equation for venous liver 
blood concentration is thus given by 


By 


Q\ 

P, 


KK- By)~ 



(5) 


The system of Eqs. (l)-( 5) represents thePBPK model for inhalation of a VOC. and 
for any given time its compartmental quantities or corresponding concentrations 
S ; or C,, are found by simultaneous numerical integration of the system. 

Exposure by routes other than inhalation is easily incorporated into this model. 
Since blood draining the stomach, small intestine, and colon passes through the he¬ 
patic portal vein, exposure to VOCs via ingestion is modeled simply by assuming a 
direct introduction of the ingested mass into the liver compartment (NAS, 1986). 
The latter introduction can be assumed to reflect a first-order infusion process, the 
approach taken in NAS (1986). or it may be modeled more simply as a constant 
infusion into the liver at a rate R ( mg/hr) by adding the constant R/0 7 , P,) to the 
right side of Fq (5). 


Pharmacokinetic Quantities of Regulatory Concern 

Of particular interest, in the context of regulatory risk assessment for certain sus¬ 
pected environmental carcinogens like TCE, is the metabolized fraction of the total 
quantity of chemical potentially available for absorption and metabolism. Under 
conditions of purely respiratory exposure ‘he corresponding quantity of interest is 
the fraction, / nir , of the maximum plausible metabolic rate, or metabolic clearance 
fraction, given a constant ambient concentration C in . For the Ramsey-Andersen 
PBPK model, this quantity may be defined as 
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/mr=.im/<ftC ra ). (6) 

which may be contrasted with the assumption, sometimes made in the context of 
carcinogen-risk assessment (Anderson ei ai, 1933). that 100% of a chemical entering 
the lungs through total respiratory ventilation (or that contained in approximately 
20 nr’/d^y for a reference 70-kg man) is absorbed and potentially available for me¬ 
tabolism. 

In the context of purely ingeslive exposure, a corresponding (hypothetical) quan¬ 
tity may be defined as the fraction. ./ mo , of the maximum plausible metabolic rate 
given a continuous rate. R. of ingestive absorption. Again by the Ramsey-Andersen 
model, this quantity would be defined as 

/mo = - UJ R- (7) 

Because ingestion is actually a discontinuous process with some average absorption 
rate R. the maximum plausible metabolic rate will be less than / m0 as defined by Eq. 

(7) to the extent that metabolic saturation occurs, allowing more of the dose to escape 
via exhalation and thus to avoid being metabolized. For VOCs. both./ mr and./^ 0 are 
necessarily less than one whenever ( finite I metabolism takes place, because any un- 
metabolized compound is always subject to pulmonary excretion. 


ANALYSIS OF THE RAMSEY-ANDERSEN PBPK MODEL 
AT STEADY STATE 

Very low-level, continuous exposure scenarios are typically of concern in the con¬ 
text of environmental regulation. The following is an analysis of how the Ramsey- 
Andcrsen PBPK model behaves under stcadv-state, respiratory exposure conditions. 
Steady-state ingestive exposure, using the constant infusion model described earlier, 
is also considered for the purpose of comparison and because the mixed exposure case 
represents a scenario of human environmental exposure which may be of particular 
regulatory concern, 

Because at steady state i? a = B, for / = 2. 3, 4. Eq. (2) reduces under steady-state 
conditions to 


5 = 


_1_ 


[<?■*, +ft(Gh-Gi)]. 


so that Eq. (I) reduces to 


D _ 0aC, n + Q\B\ 

a (aw+ft’ 


(8) 


(9) 


Also at steady state, Eq. (5), modified as described above to reflect constant ingestive 
infusion, reduces to the form 


' -Qi(g a -#,) + /?■ (10) 

A m + U\ 

so that the solution for venous liver blood concentration, given input C, n , is the qua¬ 
dratic root 
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( 8 ) 


(9) 


it ingestiye 


( 10 ) 

is the qua- 


1 1 


\ 


H x -- Y+VY i + Z. 

in which 


(Ha) 


and 


T = — [C^P h + (R — T n , a ,i)W— A m ], 

(Hb) 

Z= K m (C, n P*+ Rtt’). 

(He) 


H = 


—? + — 
Q, Qi 


(lid) 


Note that the parameters Qb . land P, for t ~ 1.4 and Q, for i = 1. 2, 3 do not 

appear in this solution. The steady-state metabolic rate is thus given by using Eqs. 
(1 I a)-(1 Id Ho evaluate 8 X in Eq. (4 >. and likewise in the expressions for the fractions 
,/ mr and of maximal metabolic rate I'orTCE given by Eqs. (6) and (7), respectively. 

Important in the context oi'envtronmenlal risk management is the limiting value 
of /, nr as C,„ — 0. that is. at very low exposure levels that might be typical of nonoccu- 
pational. purely respiratory exposure to VOCs. At such exposure levels, metabolism 
is essentially a linear function of applied dose in the Ramsey-Andersen PBPK model. 
In this discussion of the low-dose situation, the ratio (l r ,mu,/A m ) shall be represented 
by the linear metabolic clearance rale A (in liters/ hr). The limiting value referred to 
is given by 


lim / mr = 

£r-0 


+ ^(1 + 2 . 

Pb \ A Q\ 


( 12 ) 


whereas the maximally conservative assumption that infinite metabolism is ap¬ 
proached as steady-state C, n 0 yields the corresponding limiting value 


lim f* r 

c— ji 


r. + ( Qt! Q \) 
[ /J * 


(13) 


representing the physiologically determined upper bound on the fraction of respira¬ 
tory intake capable of being metabolized. This value is a function of just three param¬ 
eters, onlv one of which is influenced by the particular chemical under consideration. 
Note that /*, and /** also apply to a VOC dose received via dermal absorption, 
since such a dose would enter systemic circulation and thus be subject to pulmonary 
excretion in the same way that a respired dose would. 

The consequences of the Ramsev-Andersen model for purely respiratory exposure 
suggested by Eqs. (12) and (13) are in contrast to the assumption regarding TCE 
metabolism made in the Fernandez et al. (1977) model, discussed earlier. Recall that 
the latter model presumed a purely flow-limited metabolic clearance fraction of 
f c = 0.86 (with a corresponding metabolic rate of in the notation used to 

describe the Ramsey-Andersen model). This assumption is different from the im¬ 
plication of Eqs. (12) and ( 13 ) that metabolic clearance of a respired VOC is directly 
(but non linearly) proportional to the corresponding clearance rate parameter, K, and 
only becomes constant asA'-» -x . Indeed,, under steady-state conditions of respira¬ 
tory exposure, Eq. ( 5 ) implies that the relation 
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I 


Q\ ( B 3 B\) _ A 
~<2, ft A'+O, 


(14) 


follows from the Ramsey-Andersen PBPK. model. Note that, using Eq. (14), Eq. (12) 
may be rewritten as 


1 I- 


iQjQxj r 1 

pj* 


(15) 


which, compared with Eq. (13). clearly shows that /* r approaches its maximal 
value. /£*, only as f. approaches unity. 

The expression for / * r given in Eq. (12) may be compared to a similar limiting 
value of/ mn as a hypothetically continuous, purely ingestive dose approaches zero, It 
can be-shown that this limiting value is given by 


f 


lint .(mu 
/?—o 



<?. c >.) 



(16) 


Therefore, /*„ > ! * r . and /£„ approaches unity under the maximally conservative 
assumption that A' — x as A — 0. That is. unlike the case of purely respiratory’ 
exposure, it is physiological])’ plausible that close to 100 fl t of a low-level, continuously 
ingested dose of a VOC can be metabolized. This is true because of a “first-pass" 
effect on ingested VOC doses, whereby they are collected in blood draining the gastro¬ 
intestinal tract and transported through the hepatic portal vein directly to the liver, 
before being passed to mixed venous blood where they become subject to respiratory 
elimination. 

In the case of mixed respiratory and ingestive exposure under steady-state condi¬ 
tions, the fraction. f m , of the maximum plausible metabolic rate is the weighted aver¬ 
age 

/• _ ,/mf(3>E in T Jmt>B 

' m_ Q.Q„+R ’ 


which, as C io — 0 and R — 0, approaches the limiting value 


_ /£ r Q.O. ±/ZqR n7) 

<2.C in + R ’ 

where f m < /*. Equation (17) remains approximately true when Cj„ and R are re¬ 
placed by their corresponding time-weighted average values, 

C in = f C, n (t)dt and R = f R(t)dt. 

Jo Jfl 

provided that temporal discontinuity in Cj„(f) and R(t) does not result in a significant 
deviation from linear, nonsaturated metabolism (i.e., provided that { (l) | (Cj n (f), 

R{t), t) V K m for all l where 0 I T). To the extent that metabolic saturation 
does take place, the latter approximation will be less than the value /£ defined by 
Eq. (17). 

Note that the quantities /£,. /£*, /./ fm - and X defined above are all 
invariant with respect to body weight, so long as rates for flows and metabolic pro¬ 
cesses (e.g,, Qf Q>>. £>,, E m ai, A') are all assumed to scale to the same 0.7 power of 
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Fjc. 3. Data of Iked a t: a/. 119721 on geometric mean levels of total crichloro compounds (TTO in 
unite produced in r»l male Japanese workers exposed S hr/day, 6 days/week to different daily time* 
weighted average concentrations of trichloroethylene i TCEj in air. Numbers next to data points indicate 
the CiUmbcr of workers monitored at the corresponding concentration level. Error bars represent rtone 
geometric standard dev :aticn. Linear, least'Squares regression lire shown tor unweighted data points. 


body weight. 1 By the same reasoning, the quantitiesare also invariant 
with respect to body weight whenever the pharmacokinetic system involved is en¬ 
tirely linear, or nonsaturated. It also follows that /£„ remains invariant with body 
weight when both .i„, and R are in units of milligrams per kilogram per day. whereas 
.-i m in milligrams per kilogram per day will be smaller tor the larger of two animals 
exposed to a given very low concentration. C in . ofVOC in air by the factor (w.,/ 
iv 2 ) o i . where m ( and im are the weights of the smaller and larger animals, respec¬ 
tively. 


APPLICATION OF PBPK MODEL TO HUMAN 
DATA ON TCE METABOLISM 

Data are available on the extent to which urinary metabolites are produced in 
workers exposed to TCE in atr at concentrations of 0 to 175 ppm (lkeda el al 1972; 
Ikeda. 1977). These data, shown in Fig. 3, are derived from surveys of workplace air 


1 This approach is often taken because, in the context of pharmacokinetic modeling, such physiological 
parameters are generally assumed to vary with basal metabolic rate m proportion to body surface area ( or. 
approximately, to body weight to the 0.7 power), rather than to body weight per sc (Gehringet al.. 1978; 
Dedrick and Bischuff. 1980; Andersen el al.. 1980; Andersen. 1981; Calabrese, 1983; Ramsey and Ander¬ 
sen, 1984; U.S. EPA. 1985. !986a; NAS. 1986). For example, given an estimated tn milligrams per 
hour for an animal of weight w ,. the corresponding predicted values for animals of weight would be 
J>* ITMJ i" *f 2 /W|)“’{in mg/hr) or. on a weigh t-norm alized basis, ( L ra0 , / w,)( w.J V,, j' J {in mg/kg-hr) That 
is. a heavier animal would be expected to metabolize less per unit body weight than a lighter one in a given 
amount of time, In contrast, the value of the Michaelis constant, K m . fox a given compound, in the absence 
of data indicating otherwise, is generally assumed to be independent of body size when this constant is 
expressed as the reactant concentration (in tng/liter blood) at which the metabolic reaction rate is one- 
half its maximal value (Ramsey and Andeisen. 1984; U.S. EPA, 1986a; NAS, 1986). Since blood wetght 
vanes with body weight approximately lo the first power for animals of widely varying weights (Dednck 
and BischofT. 1989), (he independence of K m and body sice should also be expected to hold for K„, values 
expressed in terms of milligrams per kilogram body weight, If the physiological clearance parameters arc 
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and urine samples collected from 51 male workers in 10 TCE workshops in Japan. 
Workplace air concentrations of TCE were reported to be relatively constant over the 
common 8 hr/day, 6 davs/week occupational schedule of the workers studied. Urine 
samples were passed at about 1:00 p.m. in the "latter half of a week” (ikeda et al. 
1972). Corfcentrations of metabolites (TCEE, TCA, and TTC) in urine were mea¬ 
sured and plotted against time-weighted average TCE-air concentrations. The rela¬ 
tionship between the reported geometric mean TTC levels in urine (adjusted to a 
specific gravity of 1.016) and corresponding ambient TCE concentrations ranging 
from 0 to 175 ppm is well represented by a O-intercept-regression line with a slope ot 
5.0 mg TTC/liter urine/ppm TCE in air (Fig 3 }. 

A PBPK model based on the Ramsey-Andersen model was adapted, as described 
below, to model TCE pharmacokinetics in humans occupationally exposed to ambi¬ 
ent TCE at the levels and in the temporal pattern experienced by the Japanese work¬ 
ers observed by Ikeda el at. (1972). For this analysis, the Ramsey-Andersen model 
was extended to include the TCE-metabolism model used by Fernandez ct al (1977), 
discussed above and illustrated in Fig. 1. An explicit model of human urinary output 
was also included as well as a urine-sampling scheme similar to that used in the Ikeda 
et al. study. 

Reference values for the parameters listed in Table 2 and used in the PBPK model 
for TCE-cxposed workers are given in Table 3. The values for the physiological pa¬ 
rameters (>,, Qi, Q,. and I". pertaining to a reference 70-kg male worker appearing 
in Table 3 are taken from a study by Ward ct al (1988). who examined a PBPK 
model for humans exposed to perchloroethylene; these values are similar to those 
used in the N AS (1986) analysis. The values for the blood / air and fat/blood partition 
coefficients ( P b and P;) for TCE were taken from U.S. EPA (1985). as adapted from 
Sato et al. (1977). The values used for/) and P 4 are based on the liver/blood, muscle/ 
blood, and biood/air partition coefficients for rats (equal to 1.69. 0.63, and 25.82, 
respectively) reported by Sato et al. (1977), under the assumption that the corre¬ 
sponding tissue/air partition coefficients for rats and humans are equal. (Supporting 
this assumption is the fact that the fat/air partition coefficients for rats and humans 
are 661 and 674, respectively, as reported by Sato et al.). Similarly, the value used 
for P 3 is based on an average tissue/blood partition coefficient ot 1.26 for rat visceral 
tissues calculated from the data of Sato et al. Finally, the metabolic rate-parameter 
values used by Fernandez et at. (1977), and shown in Fig. I. were used to represent 
those of a reference 70-kg worker. 

It was assumed that a 70-kg male produces an average of 0.0583 liter urine per 
hour during the day, for unne adjusted to a specific gravity of 1.016 (1CRP. 1975). 
In the absence of more specific data, it was further assumed that urine samples were 
passed by the workers studied by Ikeda et al. (1972) at 1:00 p.m. on Thursdays, Fri¬ 
days, and Saturdays of each week, and that a prior urination (emptying the bladder) 
occurred 4 hr before (i.e., at 9:00 a.m.) on each collection day for all workers. 

For the purpose of modeling the occupational exposures studies by Ikeda el al. 
(1972), the reference parameter values listed in Table 3 were scaled to approximate 


assumed to be directly proportional lo the 0.7 power of body weight, it follows that kinetic rate constants 
governing uptake, metabolism, and excretion—all in uniis ol per hour, or liters per houT per liter compart¬ 
ment volume where 1 liter body volume is assumed to weigh about l kg—are i nversclv proportional to the 
0.3 power ofbody weight, 
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Par -.meter Values Usedin PBPK Model for TCE 



Reference 

Male reference 

Japanese 

Parameter 

Unit 

rat d 

human'’ 

male worker f 

W'< weight bodv) 


0.30 

70 

55.2 

ft 

liters/hr 

5.74 

353.5 

299.3 

ft 

liters/hr 

5.74 

371.6 

314.7 



22.0 

9.92 

9.92 

e,/ft. (= i 


0.75 

0.25 

0.25 

i 


0,09 

0.05 

0.04 

3 


0.47 

0.51 

0.52 

4 


0.19 

0.19 

0.19 

T,/H\/= 1 


0.041 

0.04 

0.027 

■> 


0.09 

0.20 

0.15 

3 


0.059 

0.05 

0.10 

4 


0.73 

0.6: 

0.61 

P„ : ■■■ 1 


1.3 

4.4(1 

4.40 

T 


26 

68 0 

68.0 

3 


13 

3.28 

3-28 

4 


0.5 

1.64 

1.64 

1 

mg/hr 

5.17 

— 

— 


mg/1 iter 

0.35 

- 

- 


■' Parameter values taken from NAS metabolic parameters were reported to be determined by 

gas uptake experiments (NAS. 1986). 

* Physiological parameter values taken from Ward ei at. (1988); values for partition coefficients are based 
on those listed in U.S. EPA (1985), derived front Sato et at. (1977). 

‘‘ Parameter values are scaled from those corresponding to 2 reference male human, based on the anatom¬ 
ical model of Kerr (1979) for reference Japanese adults. 

those applicable to a typical Japanese male worker. To this end, body weight was set 
at 55.2 kg and the tissue volumes V, were adjusted to the new values shown in Table 
3, based on an anatomical model for reference Japanese adults developed by Kerr 
(1979). As shown in Table 3. the values for the blood-flow fractions QilQb and Q 3 / 
Q b were changed slightly to reflect the altered tissue volumes. The reference flow rates 
0 a and (2b and that for urinary output were ail decreased by the factor ( 55.2/70) 07 . 
and the reference rate constants for metabolism (appearing in Fig. 1) were all in¬ 
creased by the factor (70/55.2)° 3 . for the reasons discussed at the end of the previous 
section. 

Workers were assumed to be exposed to TCE from 8:00 a.m. to 12:00 noon and 
1:00 to 5:00 p.m. (i.e., 8 hr/day with a 1-hr lunch break) on Monday through Satur¬ 
day of each week. To approximate conditions of dynamic equilibrium, the occupa¬ 
tional exposure scenario was run for a simulated 5-week period. Numerical integra¬ 
tion of the system of differential equations involved was done on a VAX 11/750 
computer using a variable-step Gear method (Hindmarsh, 1983). On each of the last 
3 urine collection days of the simulation, calculated urinary concentrations of TCA 
and TCEL (adjusting for the molecular weight difference between these metabolites 
and TCE) were added and the three sums were averaged to yield a predicted concen¬ 
tration of TTC in unne corresponding to any given input value for the metabolic 
clearance parameter K = (F ma ,/K m ). 
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The PBPK model described was used to estimate a unique value of K that allowed 
the model to predicta metabolic output of 5.0 mg TTC/ liter urine/ppm TCE, consis¬ 
tent with the data ofIkeda et (tl. (1972) for workers exposed to TCE in the manner 
described above. The estimate of K we obtained is 47.000 liters/hr, which, from Eq. 
(14). implies that 99.8%. i.e.. virtually all. the TCE in blood entering the liver is 
metabolized in humans so exposed. 


APPLICATION OF PBPK APPROACH TO CaNCER-RISK 

ASSESSMENT FOR ICE 

Metabolized Dose in !lunutns 

On the basis of the foregoing analysis of TCE metabolism in humans, a good esti¬ 
mate of the metabolized fraction f & r oFTCE inhaled by humans at ambient concen¬ 
tration levels of0 to j 75 ppm (and ofdcrmally absorbed TCE at similarly low expo¬ 
sure levels) is provided bv Eq. (1 3 ) to be about 72%. This value of 72% is somewhat 
less than the mean value of 91‘7 (range 87 to 93%) for the “percentage of retained 
dose metabolized" (PRDM) based on the data from three studies of TCE metabolism 
in experimentally exposed humans discussed earlier and summarized in Table 1. This 
discrepancy is explained by the fact that, by definition. t m 4 PRDM. because PRDM 
already takes into account that only a (time-dependent) fraction.,/ r . of respired air is 
cleared of TCE to constitute the "retained" TCE dose referred to by PRDM, such 
that / mr = PRDM X f t . If it is assumed that the average value of/i was approximately 
0.78 for the three studies referred to (see Table 1, footnote b). then the. product 
PRDM X f, equals about 71%. which is consistent with the estimate of 72% de¬ 
rived above. 

Again using the estimate of A' obtained above to describe the Ikeda et al. data, the 
fraction f„ 0 of ingested TCE that is metabolized by humans is given by Eq. (16) to 
be approximately 100%. Consequently, for humans subject to low levels of TCE ex¬ 
posure by both respiratory and ingestion pathways. Eq. (17 ) may be rewritten as 

24H m = 24(0.72(2 a Ci n 4- R) (18a) 

or 

24zi m = l7<2 a C,« + 2C W , (18b) 

which represent the daily dose of TCE metabolized (in mg/day). given a daily, time- 
weighted average ambient concentration C, n (in mg/litejr) and ingestion rate R (in 
mg/hr) or water concentration C* (in mg/liter = 12/f, assuming ingestion of 2 
liters/day). 


Carcinogenic Potency Estimation 

To assess potential increased lifetime cancer risk (Il.CR ) to humans from known 
rodent carcinogens in the absence of adequate epidemiological data, risk may be ex¬ 
trapolated on the basis of tumor response data from rodent bioassays (Andersen et 
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Fig. 4. Simulation of metaholism ol’triehloroethylene (TCE) in a 0.30-kg rat exposed via (a) inhalation 
and(bl ingestion, using the Ramsey-Andersen PBI’K model with the parameter values listed in Table 3. 


a!.. 1983; U-S. EPA. ! 986b of which several have been carried out in the case of 
TCE (NCI. 1976: Belief al.. 1978: Henschler et al.. 1980; NTP, 1983; Fukuda, 1983; 
Maltoni e 1 al.. 1986). However, for metabolically activated compounds like TCE, 
eachapplied, lifc'ime time-weighted average (L TWA) dose rate used in the bioassays 
must first be converted to a corresponding effective dose rate. A/, ot T CE metabolites, 
either in milligrams per kilogram per day or milligrams per square meter per day 
(using, respectively, a body-weight or surface-area approach to interspecies extrapola¬ 
tion of doses assumed to be equipment with respect to eliciting a tumorigenic 
response). Data from mass-balance studies of TCE metabolism in rodents dosed with 
radiolabeled TCE may be used to convert D to A/, provided that the test species and 
strains are the same and that either (i) the exposure regimens and the animal weights 
used in the metabolite studies and the bioassay are identical or (ii) the pharmacoki¬ 
netics involved are sufficiently invariant with respect to exposure duration and body 
weight that/).,,, and / mo are not expected to vary from one exposure scenario to another 
tor animals of differing weight. 

Simulations using the Ramsey-Andersen PBPK model, with the parameter values 
for a reference rat listed in Table 3, predict that rats (and hence mice, which weigh 
less than rats) dosed with TCE by gavage or inhalation at the levels and in the tempo¬ 
ral patterns used in the TCE cancer bioassays cited metabolize 100% of retained TCE 
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TABLE 4 

Bioassay Dose-Rfsponse Dvr*, and Corresponding Estimates 




of Carcinogenic Potency for I CE 




Study, spec, 
strain 

MaKom et « 
1986. mi* 

Study, species, 
strain 

Sex, 

weight 

[dosed 

animals) 

Daily 

experimental 
applied dose or 
Concn.. D 

I.TWA 
metabolized 

tinsel .17 
Img/kg-day) 

Tumor 

Type* Incidence' 

95% UCL potency * | 

of metabolized de w: j 

(mg .V/kg-day)" 1 ", 

BW' SA' 








V 

Swiss 

NCI 1976. 

M. 34 g 

0 mg/kg 

0 


1/20 




mice. 


1169 mg/kg 

369.6 

HCC 

26/48 

0.0025 

0.032 

_ _ _ _ 

B6C3F1 


2339 mg/kg 

739.4 


31/40 


1. 

Lifetii 






FT -IT 

0.0033 

0.042 V 







TT-LT 

0.00034 

0.0044 

on deriva 


F.39g 

0 mg/kg 

0 


0/18 



h HCC. 



869 mg/kg 

274."’ 

HCC 

4/42 

0.00073 

0.0098 , 

noma; R 



1 1 39 mg,'kg 

549 8 


11/37 



malignan 









‘ Tumi 

NTP 1983. 

M 37 g 

0 mg/kg 

0 


8/48 




mice. 


1000 mg/'kg 

56.3 

HCC 

30/50 

0.0019 

0.023 

tumor t> 

B6C3F1 





TT-tT 

0.0029 

0.036 > * 

“inciden 1 






TT-LT 

00013 

0.016 ' 

d ’Pob 



0 mg/'kg 

0 

HCC L>r 

l 1/48 



stage coe 



KXXi mg/ kg 

563 

HC\ 

38/50 

0.0029 

0.036 

time-io-l 






rr-iT 

0.0048 

0.059 L 

1979;U 






rr-LT 

0.0018 

0.022 

UCL, or 


F. '-.H 

o mfc- kg 

n 


2/41 



' BW. 



1 (XX) mg/ kg 

563 

HCC 

13/4 1 

0.00096 

0.012 

SO .C/hum. 






TT-1T 

0.0011 

0.014 <1 

'SA, s 






TT-LT 

000077 

0.0099 

SO A/foym 



0 mg/kg 

0 

HCC or 

4/41 





1000 mg/kg 

563 

HCa 

19/41 

0.0014 

0.018 







TT-IT 

0.0019 

0.024 







TT-LT 

0.00 1 1 

0.014 










within 

NTP 1983, rats. 

M. 340 g 

0 mg/kg 

0 


0/33 




F344/N 


500 mg/kg 

198 

RTC 

0/20 

0.00074 

0.0043 

repeat 



1000 mg/kg 

282 


3/16 



on TC 



0 mg//kg 

0 


0/45 


i 

180 d 



500 mg/Vg 

198 

RTC or 

2/39 

0.00065 

0.0038 



1000 mg/kg 

282 

RTA 

3/26 



respet 






TT-IT 

0.00034 

0.0020 

data c 






TT-LT 

0.00068 

0.0040 
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Bell et ai 1978, 

M, 35 g (?) 

0 ppmv-6 hr 

0 


18/99 



data t 

mice. 


100 ppmv-6 hr 

42.3 


28/95 



urate. 

B6C3FI 


300 ppmv-6 hr 

[27 

HCC 

31/100 

0.0020 

0.026 



600 ppmv-6 hr 

254 


43/97 


% 

revea 



0 ppmv-6 hr 

0 


20/99 



of fm 



100 ppmv-6 hr 

42.3 

HCC or 

35/95 



i n XeV 



300 ppmv-6 hr 

127 

HCA 

38/100 

0.0028 

0.036 

111 a Cl 



600 ppmv-6 hr 

254 


53/97 


'■ — 

TCE 

Henschlcr et at 

F, 30 g (?) 

0 ppmv.6 hr 

0 


9/29 



says i 

1980, mice, 


100 ppmv-6 hr 

33.2 

ML 

17/30 

0.0074 

0.098 


Han;NMRl 


500 ppmv-6 hr 
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18/28 


i 
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Fukuda el ai 

F,30gn 

0 ppmv-7 hr 

0 


J/49 




1983, mice. 


50 ppmv-7 hr 

25.8 


3/50 




ICR 


150 ppmv-7 hr 

77.4 

LA 

8/50 

0.0014 

0.019 

Bt 
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232 


7/46 
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potency* 
lized dose 
W) 

g-day) -1 


Study, species, 
strain 

Maltom ft al 
1986. mice. 
Swiss 


Sex. 

wei&hi 

(dosed 

animals) 

Daily 

experimental 
applied dose or 
Ooncn.. D 

LTWA 
metabolized 
dose". 1/ 

(mg/kg-das > 

Tumor 

Type* incidence 4 ' 

M. 41 * 

0 ppmv-7 hr 

0 

4/90 


100 ppmv-7 hr 

J5.3 

2/90 


300 ppmv-7 hr 

;C6 

MH 8/90 


600 ppmv-7 hr 

212 

13/90 


95^ UCL potency* 
nf metabolized dose 
$;<-«) 

(mg Vf/kg-day)' 


0.00082 0.0093 


* Lifcri rite, time-weighted-average metabolized dose. .if. in ir.g/kg-day. See Bogen et til, (1988) lor details 
ort denvation as a function of P. 

* HCC. hepatocellular carcinoma; HCA, hepatocellular adenoma; RTC. renal tubular-cell adenocarci¬ 
noma; RTA, renal tubular-cell utlenoma, ML. malignant lymphoma; LA, lung adenocarcinoma: MH, 
malignant hepatoma. 

: Tumor-ircidenec denominator excludes animals dying before the occurrence of the first corresponding 
lumor type observed in the NCI 1 1976) and NTP (1983) studies TT-1T. time-to-mmor data using an 
"incidental-tumor" model; TT-LT. time-to-tumordata using a ‘lethai-tumor" model. 

* ‘‘Potency'' here means the low-dose dose-response slope expressed b> an upper-bound linear multi¬ 
stage coefficient such that at very low doses, risk - t potency x dascl. according to a multistage (or, \ritb 
Urne-io-tumor data as input, a time-dependent multistage) nsk-predicnon model (Crump and Watson, 
1979; U S. ERA. 1980; Anderson el al.. 1983; Howe and Crump, 1983: Crump and Howe, 1984). 95% 
UCL. one-tailed 95't upper confidence limit. 

r BW. body weight interspecies dose-cxtrapolatmn method, equivalent doses assumed to he in nryg/kg, 

so • 

'S A. surface area jnterspecies dose-extrapolation method: equivalent doses assumed to be in mg/kg‘T 
so A/ hum „ - -'/'m.mu((animal weight)/70 kg|‘ 


within 24 hr (Fig, 4), Thus, no buildup of residual TCE within body tissues upon 
repeated daily dosing is expected. This prediction is supported by experimental data 
on TCE metabolism in B6C3F1 mice given daily gavage doses of TCE tor 1, 10, and 
180 days (Green and Prout, 1985). Recalling that /* 0 and f%, are invariant with 
respect to body weight, it follows, under conditions of linear metabolic kinetics, that 
data on TCE metabolism from single-exposure mass-balance studies may be extrapo¬ 
lated directly to the rodent-bioassay context (correcting for nonexposure days). The 
data of Prout el al. (1985) dearly demonstrate that TCE metabolism in mice is unsat¬ 
urated even at very high levels of TCE exposure. Furthermore, PBFK. simulations 
revealed that small changes in rodent body weight do not influence predicted values 
of / mo and / mr . even at high levels of TCE exposure involving metabolic saturation 
in rats. Thus, experimental metabolism data may be used to predict the extent of 
TCE metabolism in rats and mice under the exposure conditions of the TCE bioas¬ 
says cited. 

Prediction of Human Cancer Risk 

Based on the rodent-bioassay data for TCE carcinogenicity cited above and on 
mass-balance studies of TCE metabolism in B6C3FI mice and Osbome-Mendel rats 
exposed via single gavage or inhalation administrations of TCE (Prout et al., 1985; 
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Stott et a/, 1982). upper-bound carcinogenic potencies forTCE in rodents were cal¬ 
culated (Bogen el a!.. 1988) and are listed in Table 4. Each of the 5D alternative 
potency values shown is a function of \f and tumor-incidence data for rodents in a 
given bioassav. The values range from0.00034 to0.098 (mg/kg-dayr ! (>290-fold 
range). These values may be used directly with estimated values of M for humans in 
order to derive estimates of (upper-bound) ILCR associated with exposure to TCE. 
For example, if a LTWA exposure to 1 fig! m 1 ( 5.38 ppb) TCE in air and 1 ug/litcr 
TCE in water is assumed, along with LTWA alveolar ventilation and water consump¬ 
tion rates of 5.05 lilers/kg-hr (from Table 3) and 2 !iters/70 kg-day, respectively, 
then Eq. (18b I and the potencies from Table 4 yield a corresponding predicted ILCR 
between 3.9 X 10 8 and l.J X l0'\ 


DISCUSSION 

Application of a PBPK approach to the estimation of VOC metabolism in humans 
may provide simple, convenient relationships between applied and toxieologicallv 
effective VOC doses under low-level exposure conditions pertinent to environmental 
regulation. Parameterization*; of PBPK. model of TCE metabolism in humans, using 
human data on TCE metabolism alter chronic exposure to between 0 and 175 ppm 
TCE in workplace air. predict that approximately 10()5r of TCE entering the liver, or 
about 72 * ol alveolarJy respired TCE. is metabolized. This prediction in accordance 
with experimental data on i CE metabolism obtained using humans exposed to TCE 
for short periods. The linear metabolic kinetics observed forTCE metabolism at oc¬ 
cupational exposures of up to 1 75 ppm TCE differ from those predicted for tetrachlo- 
roethylene using a similar PBPK. approach: the latter kinetics have been shown to 
involve saturation at occupational exposure levels between 10 and 100 ppm in air 
(Bogen and McKotie, 1987 i. Nevertheless, linear pharmacokinetics are expected for 
both compounds, and probably for many other halogenated VOCs as well, at the very' 
low exposure levels of environmental regulatory' concern. This fact often allows a 
simplified PBPK approach, of the type outlined here, to be adopted in addressing 
problems of regulatory' toxicology such as environmental carcinogen-risk assessment 
for VOCs. 
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